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Key  Points: 

•  Observations  showed  complex  flow 
over  the  East  Flower  Garden  Bank  in 
June  201 1 

•  Flow  was  blocked  and  diverged 
upstream  of  the  bank 

•  A  wake  with  elevated  turbulence  and 
mixing  developed  downstream  of 
the  bank 
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Abstract  In  June  2011,  a  9  day  oceanographic  survey  was  conducted  over  the  East  Flower  Garden  Bank, 
a  coral  reef,  located  on  the  outer  shelf  in  the  northwestern  Gulf  of  Mexico.  Current  temperature,  conductiv¬ 
ity,  and  microstructure  measurements  were  collected  to  characterize  flow  evolution,  turbulence,  and  mixing 
over  the  bank.  During  the  experiment,  the  flow  was  highly  stratified,  subcritical  (Froude  number  below  0.4), 
hydrostatic,  and  nonlinear  with  rotational  effects  being  important.  Observations  showed  that  flow  structure, 
turbulence,  and  mixing  were  highly  dependent  on  the  direction  and  strength  of  the  current;  thus,  they  var¬ 
ied  spatially  and  temporarily.  Responses  resulting  from  interactions  between  the  free-stream  flow  and  the 
obstacle  were  significantly  different  on  the  upstream  and  downstream  sides  of  the  bank.  Blocking  and 
diverging  of  the  flow  just  below  the  bank  height  was  observed  on  the  upstream  side.  On  the  downstream 
side,  a  wake  with  imbedded  vortices  developed.  Moreover,  turbulence  was  amplified  over  the  bank  top  and 
on  its  downstream  side.  Turbulent  dissipation  rates  were  as  high  as  10  6  W  kg  1  and  resulted  in  measured 
rates  of  energy  dissipation  and  mixing  by  turbulence  per  unit  width  as  high  as  40  W  m  \  Mixing  on  the 
downstream  side  was  elevated  with  eddy  diffusi vities  reaching  10  3  m2  s  ',  well  above  a  typical  value  of 
10  5  m2  s  1  commonly  found  in  the  ocean  thermocline  and  over  shelves  with  flat  topography.  On  the 
upstream  side,  estimated  eddy  diffusivities  were  close  to  that  for  the  ocean  thermocline,  i.e.,  they  were  gen¬ 
erally  less  than  0.5X10  4m2s  \ 


1.  Introduction 

Substantial  research  efforts  in  past  decades  have  been  devoted  to  studying  flows  over  rough  bottom  topog¬ 
raphy  in  the  ocean.  They  have  led  to  the  conclusion  that  mixing  is  much  greater  there  than  mixing  over  a 
gradually  varying  seabed  or  in  the  thermocline  in  the  open  ocean  [e.g.,  Gregg,  1 987;  Dewey  and  Crawford, 

1 988;  Ledwell  et  al.,  1 993;  Lueck  and  Madge,  1997].  Significant  bathymetric  features  such  as  sills,  canyons, 
banks,  seamounts,  ridges,  underwater  volcanos,  and  mountain  ranges  are  ubiquitous  in  the  ocean.  They 
impact  flow  patterns  locally,  and  they  may  even  influence  large-scale  circulation  [e.g,  Kunze  and  Toole, 

1 997;  Polzin  et  al.,  1 997;  Nash  and  Mourn,  2001 ;  Rudnick  et  al ,  2003;  McCabe  et  al.,  2006], 

On  continental  shelves,  small  banks  are  fairly  common  features.  Their  spatial  dimensions  can  vary  greatly 
from  a  few  hundred  meters  to  several  kilometers;  however,  they  are  often  too  small  to  be  included  in 
numerical  models  developed  for  the  coastal  waters.  As  a  result,  flow,  stratification,  mixing,  and  biochemical 
processes  are  not  well  resolved  by  numerical  simulations  near  such  bathymetric  structures.  In  situ  observa¬ 
tions  have  shown  [e.g,  Loderetal.,  1992;  Mourn  and  Nash,  2000;  Nash  and  Mourn,  2001;  Dewey  et  al,  2005] 
that  underwater  banks  significantly  alter  ambient  flow,  considerably  increase  density  gradients,  modify  ver¬ 
tical  current  shear,  and  enhance  turbulence.  Stratified  flow  approaching  an  obstacle  often  accelerates,  and 
to  conserve  momentum,  mechanical  energy  is  released  that  can  be  dissipated  locally  or  radiated  away 
through  internal  waves.  As  a  result  lee  waves,  hydraulic  jumps,  and  eddies  are  very  common  features  near 
banks.  Amplified  turbulence  over  a  bank  usually  accounts  for  a  large  fraction  of  the  entire  shelf  mixing, 
exerts  a  larger  bottom  stress,  and  dissipates  a  significant  amount  of  energy  in  a  shelf  region  [e.g.  Mourn 
and  Nash,  2000]. 

Stratified  flows  over  an  obstacle  are  often  categorized  as  subcritical,  critical,  and  supercritical  in  terms 
of  the  Froude  number  that  is  defined  as  a  ratio  of  the  flow  speed  to  the  buoyancy  frequency  and  the 
total  depth.  Evolution  of  flows  and  their  major  features  vary  dramatically  as  a  function  of  the  Froude 
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number.  Past  observational 
and  modeling  efforts  in  the 
ocean  have  often  been 
focused  on  studying  ener¬ 
getic  flows  with  Froude  num¬ 
bers  approaching  1  or  larger, 
i.e„  critical  and  supercritical 
flows  over  rough  topography 
with  energy  being  released 
through,  for  instance,  a 
hydraulicjump  [e.g.,  Wesson 
and  Gregg,  1 994;  Farmer  and 
Armi,  1999;  Hogg  etal.,  2001; 
Mourn  and  Nash,  2000].  Much 
of  our  understanding  of  flow 
dynamics  for  very  low  Froude 
numbers  (much  less  than  1) 
comes  mostly  from  the 
atmospheric  research  of  sub- 
critical  airflows  over  a  moun¬ 
tain  [e.g.,  Hunt  and  Snyder, 

1 980;  Smolarkiewicz  and  Rotunno,  1 989;  Smith  and  Grands,  1 993;  Vosper  et  a /.,  1 999;  Dalziel  et  at.,  2011], 
These  airflows  are  usually  characterized  by  blocking,  splitting,  and  moving  around  on  the  upstream  side 
of  the  obstacle.  On  its  downstream  side,  a  wake  develops  where  enhanced  turbulence  and  eddy  shed¬ 
ding  are  commonly  observed. 
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Figure  1.  Map  of  the  East  Flower  Garden  Bank  with  depth  contours  (gray  lines)  every  5  m; 
also  shown  are  locations  of  the  long  term  Barny  (Ml  M4;  blue  circles)  and  string  (SI  S4; 
cyan  diamonds)  moorings,  short  term  Bamy  (M4a,  M5a,  M6  M8;  red  circles)  and  string  (S5 
and  S8;  black  diamonds)  moorings,  VMP  drops  (blue  dots),  selected  ScanFish  tows  (black 
lines),  and  meteorological  buoy  42047. 


In  this  paper,  recent  current  and  hydrographic  observations  are  analyzed  to  describe  the  evolution  and 
main  features  of  the  subcritical  highly  stratified  flow  over  an  isolated  bank  located  on  the  outer  north¬ 
western  Louisiana-Texas  shelf  in  the  Gulf  of  Mexico.  Moreover,  collected  microstructure  measurements 
allow  an  assessment  of  turbulence  and  mixing  in  this  region.  The  paper  is  organized  as  follows.  In 
section  2,  the  bank,  field  experiment,  and  observations  are  described.  Background  meteorological  and 
oceanographic  conditions  are  discussed  in  section  3.  Snapshots  of  flow  patterns  and  stratification  are 
presented  in  section  4.  Turbulence  and  mixing  conditions  over  the  bank  are  characterized  in  section  5. 

In  section  6,  a  discussion  of  the  results  in  terms  of  nondimensional  parameters  is  presented.  Finally,  find¬ 
ings  are  summarized  in  section  7. 


2.  East  Flower  Garden  Bank  Experiment 

The  field  experiment  was  a  part  of  the  MORT  ('Mixing  over  Rough  Topography')  and  COB  ("Currents  over 
Bank")  projects  funded  by  the  Naval  Research  Laboratory  (NRL)  and  the  Bureau  of  Ocean  Energy  Manage¬ 
ment  (BOEM),  respectively.  It  took  place  from  31  May  201 1  to  1 3  June  2011  over  and  around  the  East  Flower 
Garden  Bank  (EFGB).  The  EFGB  is  located  on  the  shelf  edge  in  the  northwestern  Gulf  of  Mexico  about 
1 90  km  southeast  of  Galveston,  TX  (Figure  1).  The  bank  is  a  thriving  coral  reef.  It  is  a  part  of  the  Flower  Gar¬ 
den  Banks  National  Marine  Sanctuary  managed  by  NOAA,  and  is  one  of  many  banks  formed  by  salt  domes 
on  the  Louisiana-Texas  shelf.  The  EFGB  is  asymmetric  with  rather  steep  slopes  on  its  southern  and  eastern 
sides.  It  is  about  1 0  km  long  and  6  km  wide  with  major  and  minor  axes  approximately  aligned  in  north- 
south  and  east-west  directions,  respectively.  The  bank  top  rises  from  water  depths  of  100-130  m  to  about 
18  m  below  the  sea  surface.  The  ocean  bottom  area  covered  by  the  bank  (within  a  100  m  isobath)  is  about 
40.76  km2. 

Two  research  vessels,  R/V  Pelican  and  R/V  Manta,  occupied  waters  over  the  bank  to  measure  currents, 
hydrography,  and  turbulent  mixing.  Almost  all  data  collected  by  instruments  deployed  from  or  towed  by 
the  ships  were  taken  between  4  June  201 1  and  12  June  201 1 .  Both  vessels  were  equipped  with  300  kHz 
acoustic  Doppler  current  profilers  (ADCPs)  to  collect  current  velocity  profiles  with  vertical  resolution  of  4  m. 
Velocity  observations  were  processed  using  the  University  of  Hawaii's  CODAS  system  [Firing  etal.,  1995]  to 
correct  for  ship  motion  and  averaged  over  4  min  intervals.  A  towed  undulating  vehicle,  ScanFish  MK II, 
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manufactured  by  EIVA  of  Denmark  and  equipped  with  a  Sea-Bird  Electronics  SBE49  FastCAT  was  used  to 
sample  hydrography  from  the  R/V  Pelican.  Conductivity,  temperature,  and  microstructure  observations 
were  collected  with  a  Vertical  Microstructure  Profiler  (VMP;  the  profiler  was  manufactured  by  the  Rockland 
Oceanographic  Services  Inc.,  Canada  and  is  similar  to  a  profiler  described  by  Wolk  et  ai  [2002])  from  the  R/V 
Manta.  The  VMP  profiler  carried  a  thermistor  (FP07),  a  microconductivity  sensor  (Sea-Bird  Electronics 
SBE7-38),  two  velocity  shear  probes,  a  high-resolution  pressure  sensor,  accelerometers,  and  fine-scale  exter¬ 
nal  conductivity  (Sea-Bird  Electronics  SBE4)  and  temperature  (Sea-Bird  Electronics  SBE3)  sensors.  VMP  pro¬ 
files  extended  from  a  few  meters  below  the  surface  to  near  the  bottom.  Moreover,  turbulent  kinetic  energy 
(TKE)  dissipation  rates  (b)  were  estimated  from  small-scale  velocity  shear  following  a  formula  for  isotropic 

turbulence  [Oakey,  1982]:  s  -y  v(|£)2,  where  v  is  the  kinematic  viscosity,  and  (^-)2  is  the  variance  of  the 
current  shear.  The  velocity  shear  variance  was  calculated  by  integrating  the  spectrum  for  the  wave  number 
range  from  2  to  30  cpm  (cycles  per  meter),  while  the  spectrum  was  estimated  in  bins  of  1 024  shear-probe 
data  records,  with  a  512-record  overlap.  Selected  ScanFish  sections  and  locations  of  VMP  drops  are  dis¬ 
played  in  Figure  1. 

Additionally,  two  dye  releases  were  conducted  from  the  R/V  Pelican  to  visualize  flow  patterns  over  and 
around  the  EFGB.  The  first  release  was  on  the  western  side  of  the  bank  around  21:00  UTCon  6  June  201 1. 
The  dye  was  released  near  the  top  of  the  thermocline  (~30  m).  The  second  dye  patch  was  also  released 
near  the  depth  of  ~  30  m  but  on  the  eastern  side  of  the  bank  at  about  15:00  LfTC  on  8  June  2011.  The  dye 
used  in  the  experiment  was  a  water  soluble  fluorescein  dye  (Keyacid  Uranine  K  Liquid)  referred  as  "Uranine* 
(www.dyes.com).  The  density  of  the  liquid  Uranine  (40%  concentrate)  was  1 1 1 1  kg  m  3.  The  dye  concentra¬ 
tion  was  reduced  to  12%  to  obtain  a  target  density  close  to  1023.79  (±0.01)  kg  m  3  by  adding  662.4  L  of 
freshwater  and  45.4  L  of  isopropyl  alcohol.  The  dye  was  released  to  the  sea  through  a  diffuser  (of  about  1  m 
length)  attached  to  the  end  of  a  hose.  The  dye  volume  of  ~ 1 022. 1  L  was  pumped  at  37.9  L  min  1  while  the 
R/V  Pelican  streamed  at  a  speed  of  about  1ms  '  in  the  north-south  direction.  The  depressor  undulated 
approximately  ±2  m  from  the  30  m  depth.  We  planned  for  a  dye  patch  of  1  m  wide,  2  m  high,  and  1 .1  km 
long  with  the  initial  dye  concentration  of  14  mg  L  1  (or  14  ppm).  The  dye  was  sampled  by  fluorometers 
attached  to  the  ScanFish  and  a  REMUS  Autonomous  Underwater  Vehicle  (AUV)  [e.g.,  Moline  et  ai,  2005],  and 
by  a  WETLabs  ac-9  absorption  and  attenuation  meter  (a  490  nm  absorption  channel)  installed  on  the  Scan¬ 
Fish.  Immediately  after  dye  release,  the  Scanfish  was  deployed  in  a  lawn  mower  pattern  while  undulating 
between  10  and  40  m  water  depths  to  capture  the  dye  distribution.  Several  hours  later,  the  REMUS  AUV 
was  deployed  for  a  duration  of  10  h.  The  ScanFish  fluorometer  was  calibrated  to  measure  a  maximum  con¬ 
centration  up  to  400  ppb  with  0.05  ppb  resolution  so  that  it  was  initially  out  of  scale  since  the  dye  concen¬ 
tration  was  greater  than  400  ppb;  hence,  we  used  the  observed  linear  correlation  between  the  uranine  and 
chlorophyll  to  recapture  the  high  values  of  the  missing  dye  concentration. 

Ten  Barny  (trawl-resistant  pods)  and  six  string  moorings  were  in  place  around  and  on  the  EFGB  during  the 
field  experiment  (Figure  1).  Five  Barny  (M1-M5)  and  four  string  (S1-S4)  moorings  were  already  deployed  for 
the  COB  project  and  remained  in  the  water  for  approximately  1  year  [Teague  et  ai,  201 3].  Another  five  Barny 
(M4a,  M5a,  M6,  M7,  and  M8)  and  two  string  (S5  and  S8)  moorings  were  deployed  just  for  2  weeks  spedfically 
for  the  experiment  as  a  part  of  the  MORT  project.  All  Bamy  moorings  contained  Teledyne  RD  ADCPs  and 
either  wave-tide  gauges  (Sea-Bird  Electronics  SBE26)  or  high-resolution  pressure  sensors  called  Ppods  [Mourn 
and  Nash,  2008].  String  moorings  were  equipped  only  with  temperature/conductivity/pressure  (TCP)  sensors, 
except  for  S8  which  also  had  a  downward-looking  ADCP.  The  number  of  TCP  sensors  on  the  string  moorings 
varied  between  5  and  12,  and  they  were  either  Sea-Bird  Electronics  MicroCAT  (SBE37)  or  In  Situ  Aqua  Troll 
instruments.  Detailed  information  such  as  deployment  times  and  depths,  sampling  rates,  ADCP  frequency, 
and  bin  depth  resolutions  for  all  moorings  deployed  for  the  COB  program  is  in  Teague  et  ai  [201 3],  Similar 
details  for  the  MORT  moorings  are  given  here  in  Table  1.  ADCPs  for  the  MORT  project  were  300  kHz  Teledyne 
RDI  units,  except  for  moorings  M5a  and  S8,  which  were  furnished  with  600  kHz  Teledyne  RDI  current  profilers. 

3.  Atmospheric  and  Oceanographic  Settings 

3.1.  Atmospheric  Conditions 

Atmospheric  pressure  and  winds  from  nearby  NDBC  buoy  42047  are  shown  in  Figure  2.  Atmospheric  pres¬ 
sure  (Figure  2a)  began  increasing  on  31  May  201 1  and  remained  high  with  a  maximum  of  10225  hPa 
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Table  1.  Mooring  Summary 


Water 

Start  Date 

Sampling  Rate  (sV 

Depth  Range  or 

Mooring 

Latitude  (N) 

Longitude  (W) 

Depth  (m) 

End  Date 

Bin  Size  (m)a 

Instrument  Depth  (m)b 

Bamy  Moorings 

M4 a 

27” 56.942' 

93” 34 563' 

106 

2  Jun  201 1 

1 3  Jurr  201 1 

60/2 

10  102 

MSa 

27” 55 .983' 

93° 36.1 46' 

46 

1  Jun  201 1 

12 Jun  2011 

60/1 

3.7  44.7 

M6 

27” 56.187' 

93° 38633' 

101 

31  May  201 1 

12  Jun  2011 

60/4 

12  96 

M7 

27”  52.830' 

93”  35  806' 

120 

31  May  201 1 

12 Jun  2011 

120/4 

8  112 

M8 

27”55.983' 

93”34j025' 

112 

31  May  201 1 

13  Jun  2011 

120/4 

8  104 

TCP  Moorings 
SS 

MC1e 

27”  56.0 13' 

93”36j059' 

47 

3  Jun  201 1 

30 

5.6 

MC2 

1 1  Jun  201 1 

30 

17.3 

MC3 

30 

26.0 

MC4 

30 

36.3 

MC5 

30 

45.9 

S8 

MCI 

27”  55.998' 

93”33953' 

112 

3  Jun  201 1 

30 

6.4 

MC2 

1 1  Jun  201 1 

30 

16.7 

MO 

30 

38.5 

MC4 

30 

49.7 

MC5 

30 

60.7 

MC6 

30 

71.7 

MC7 

30 

7a7 

MC8 

30 

87.7 

MC9 

30 

96.7 

MC10 

30 

110.0 

ADCPd 

60/1 

75.7  111.7 

73j6e 

“Sampling  Rate  (s)  sampling  rates  fot  all  sensors  used  for  the  MORT  project;  Bin  Size  (m)  vertical  resolution  for  ADCP  instalments 
only. 

"Depth  Range  depth  ranges  over  which  ADCP  data  were  collected  at  Barny  moorings  and  the  downward  looking  ADCP  on  S8; 
Instrument  Depth  depths  for  MicroCATS. 
cMC  Sea  Bird  MicroCAT  SBE37  instruments. 
dADCP  the  downward  looking  RDI  600  kHz  ADCP  on  S8. 

*A  depth  of  the  downward  looking  ADCP  on  S8. 


during  the  first  part  of  the  experiment.  After  5  July  201 1,  it  dropped  below  1 016  hPa.  A  few  days  before  the 
experiment  began,  the  winds  were  fairly  strong  reaching  10  ms  1  (Figure  2b)  and  blew  from  the  southern 
quadrant,  i.e,  from  southeast,  south,  and/or  southwest  directions.  Winds  from  the  southern  quadrant  are 
typical  for  summer  months  in  the  northwestern  Gulf  of  Mexico  [Gutierrez  de  Velasco  and  Winant,  1996].  Dur¬ 
ing  the  sampling  time  period,  winds  weakened  and  dropped  below  6ms  '.Wind  directions  were  variable 
directions  until  8  July  201 1 .  Then  they  switched  again  to  a  summer  pattern,  i.e.,  to  southeasterlies. 

3.2.  Oceanographic  Conditions 

Mean  flow  (frequency  less  than  0.6  cycles  per  day  (cpd))  had  two  phases  during  the  experiment  (Figure  3). 
Figure  3  shows  40  h  low-pass  current  data.  Until  mid-day  on  8  June  2011,  mean  currents  were  moderate 
(40  cm  s  1  or  less)  and  generally  directed  eastward.  Afterward,  the  mean  currents  weakened  (30  cm  s  1  or 
less)  and  reversed,  flowing  toward  the  west.  Note  that  mean  flow  was  less  impacted  by  the  bank  at  the  Ml 
mooring  than  at  the  M4  mooring.  M4  was  deployed  closer  to  the  EFGB  and  the  currents  just  below  the 
bank  height  were  very  weak,  especially  evident  for  the  eastward  flow  phase  (Figure  3).  Higher  frequency 
fluctuations  including  tides,  inertial,  and  superinertial  oscillations  are  usually  superimposed  on  the  mean 
flow.  Barotropic  tidal  currents  are  rather  weak  in  the  northwestern  Gulf  of  Mexico  [DiMarco  and  Reid,  1998]. 
Teague  et  al.  [201 3]  estimated  amplitudes  of  major  tidal  constituents  M2,  K,,  and  O,  to  be  between  3  and 
5  cm  s  1  in  the  EFGB  region.  Near-inertial  oscillations  have  a  period  of  25.58  h  (close  to  the  O,  period),  and 
they  are  generally  more  energetic  than  barotropic  tides  in  the  northwestern  Gulf  of  Mexico  [ Chen  et  al., 

1 996].  Near  the  EFGB,  the  inertial  currents  were  about  1 5  cm  s  1  or  larger  [Teague  et  al,  201 4],  During  the 
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experiment,  total  measured  currents  com¬ 
prised  the  mean  flow,  tides,  near-inertial 
currents,  and  higher  frequency  motions 
(frequency  higher  than  2.4  cpd)  and  are 
displayed  in  Figure  4.  Tidal  currents  showed 
a  spring-neap  cycle,  while  the  inertial  cur¬ 
rents  were  likely  produced  by  rotating 
winds  at  the  beginning  of  June.  In  the 
upper  40  m,  the  currents  were  as  high  as 
50  and  25  cm  s  ',  while  40  m  depth/time 
averages  were  about  20  and  1 0  cm  s  1 
with  standard  errors  of  1  and  0.5  cm  s  1  for 
the  eastward  and  westward  mean  flows, 
respectively.  All  velocity  estimates  were 
computed  from  current  observations  col¬ 
lected  at  Ml,  which  seem  to  be  the  least 
impacted  by  a  presence  of  the  EFGB.  Only 
3.3%  of  the  EFGB  area,  located  just  in  its 
southern  part  rises  above  40  m  of  the  water 
depth  so  that  the  flow  above  40  m  should 
have  been  influenced  the  least  by  the  bank. 
This  assumption  should  be  especially  valid 
at  the  Ml  mooring  located  approximately 
7  km  north  of  the  highest  part  of  the  bank. 
Therefore,  40  m  depth  current  veloaty 
means  calculated  from  Ml  observations 
should  tie  rather  good  estimates  of  the 
free-stream  flow  velocity  and  were  used  in 
further  analyses  and  discussions. 

Time  series  of  temperature,  salinity,  and 
density  from  S8  are  plotted  in  Figure  5. 
These  observations  encompassed  almost  the  entire  time  period  of  the  field  experiment  and  showed  that 
there  was  a  very  shallow  mixed  layer  during  this  time  as  also  indicated  by  Teague  et  al.  [2013].  Salinity  at  S8 
varied  little  vertically  with  a  range  between  35 .20  and  36.66  psu,  while  temperature  showed  more  variability 
ranging  from  1 8.10°C  near  the  bottom  to  28.641>C  at  about  6  m  below  the  sea  surface.  Similar  variability  and 
ranges  were  found  at  the  other  TCP  moorings.  Such  distributions  of  temperature  and  salinity  around  the  bank 
generated  a  mean  density  that  nearly  linearly  increased  with  depth  as  indicated  in  Figure  6a.  Displayed  mean 
densities  were  calculated  from  available  TCP  temperatures  and  salinities  that  were  averaged  over  time  sepa¬ 
rately  for  each  mean  flow  phase,  i.e.,  between  00.00  UTC4  June  and  12.00  UTC8  June  201 1  for  the  mean  east¬ 
ward  flow  and  between  1 200  UTC  8  June  and  23.00  UTC  1 2  June  201 1  for  the  mean  westward  flow.  The 
density  profiles  also  implied  that  there  was  almost  no  difference  in  the  mean  stratification  between  the  mean 
eastward  and  westward  flows  over  the  EFGB.  Moreover,  the  mean  density  was  used  to  estimate  an  average 
buoyancy  frequency  (A/,  the  squared  buoyancy  frequency  is  defined  as:  N2  (g/po)(dp/dz),  where  g  is  the 
gravitational  acceleration,  p  is  the  density,  z  is  the  depth  in  meters,  and  pQ  is  the  reference  density)  shown  in 
Figure  6b.  The  N  profiles  displayed  some  variations  with  depth  and  among  TCP  mooring  locations,  but  very 
limited  variability  between  the  mean  flow  phases.  The  mean  N  values  were  generally  less  than  0.03  s  1  with 
the  depth  average  of  0.018  s  '  and  standard  error  of  0.001  s  \  Equivalent  values  for  the  upper  40  m  are 
0.016  ±  0.001  s  1  and  0.018  ±  0.002  s  1  for  the  eastward  and  westward  mean  flows,  respectively.  These  two 
40  m  N  means  were  used  in  further  analyses  as  representative  numbers  of  the  ambient  stratification. 


jb)  Winds  (m/s) 


05/25  05/26  05/27  05/28  05/29  05/30  05/31  06/01  06/02  06/03 


06/04  06/05  06/06  06/07  06/08  06/09  06/10  06/11  06/12  06/13 

i 

5  m/s 


Figure  2.  (a)  Atmospheric  pressure  (hPa)  and  (b)  winds  (m  s '  ’)  a  few  days 
before  and  during  the  field  experiment. 


4.  Flow  Patterns  and  Stratification  Over  the  East  Flower  Garden  Bank 

Figures  7  and  8  show  current  observations  from  three  transects  (PT4,  PT5,  and  T3)  recorded  by  shipboard 
ADCPs  between  12:15  UTC  on  6  June  and  230  UTC  on  7  June  201 1.  Figure  8  also  displays  hourly  current 
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Figure  3.  Flow  at  (a  and  b)  Ml  and  (c  and  d)  M4:  east/west  (Figures  3a  and  3c)  and  north/south  (Figures  3b  and  3d)  current  components  (cm  s  ');  40  h  low  passed  data;  positive  values 
are  toward  east  and  north. 


data  at  four  depths  from  all  moorings  as  well  as  means  estimated  from  those  data  averaged  over  the  dura¬ 
tion  of  the  survey  time  period.  Moreover,  Figure  7  shows  density  along  the  three  transects.  During  these 
passes,  the  mean  flow  was  eastward.  The  currents  were  generally  directed  toward  the  east  above  the  bank 
height  on  both  upstream  and  downstream  (lee)  EFGB  sides  (Figures  7  and  8a).  Below  the  crest  on  the 
upstream  side,  the  flow  decelerated  and  diverged  horizontally  around  the  bank  as  indicated  by  current 
data  shown  in  Figures  8b-8d.  The  flow  was  more  complex  on  the  downstream  side  of  the  bank.  The  north/ 
south  velocity  structure  for  the  PT4  and  PT5  sections  suggested  a  presence  of  a  vortex  on  the  lee  side  of 
the  bank  (Figures  7b  and  7d).  They  also  indicated  that  this  vortex  with  speeds  as  high  as  20  cm  s  1  was 
cyclonic  with  a  horizontal  scale  of  about  4  km  and  a  finite  thickness  of  about  20  m.  Current  data  along  PT4 
and  PT5  transects  were  taken  within  6  h  and  were  the  most  applicable  to  be  used  for  computing  a  vertical 
component  of  the  relative  vorticity  (C  |  ^  where  u  and  v  are  horizontal  east/west  and  north/south 
velocity  components,  respectively).  Figure  9c  shows  layer- averaged  £  (H  <  40  m,  40  m  <  H<  60  m,  and 
H  >60  m,  where  H  is  the  depth).  The  vorticity  on  the  upstream  side  was  negligible  and  generally  close  to 
the  noise  level  of  the  computations.  However,  a  layer  of  positive  vorticity  exceeding  10  4  s  1  was  found  at 
depths  between  40  and  60  m  on  the  lee  side  of  the  EFGB,  where  £  was  about  a  factor  of  two  larger  than  the 
planetary  vorticity,  f  6.84X  10  5s  ’.  The  velocity  field  close  to  the  bank  was  not  well  resolved  by  the 
shipboard  ADCP  measurements,  and  therefore,  estimates  of  negative  vorticity  found  close  to  the  bank  on 
its  lee  side  were  not  very  reliable.  Positive  temperature  (T7)  and  negative  potential  density  (tre')  anomalies 
(Figures  9a  and  9b)  estimated  from  observations  collected  along  the  PT4  transect  were  prevalent  between 
40  and  60  m  and  indicated  deepening  of  the  isothermals  and  isopycnals  on  the  lee  side  of  the  bank.  These 
variations  were  partially  related  to  lee-side  vortex  as  well  as  to  internal  waves  and/or  mixing.  Anomalies  of 
F  and  <t,/  at  a  given  depth  along  PT4  were  computed  by  subtracting  transect-averaged  means.  Similar  tem¬ 
perature  and  density  variations  on  the  lee  side  were  observed  along  PT5  (data  not  shown).  Typically,  sub- 
mesoscale  cyclonic  eddies  in  the  ocean  are  associated  with  negative  temperature  anomalies,  contrary  to 
the  positive  temperature  variations  observed  here.  Observations  along  the  T3  transect  indicated  highly  vari¬ 
able  currents  on  the  lee  side  suggesting  a  plausible  eddy-like  feature  below  40  m  (Figures  7e,  7f,  and  8d). 
Furthermore,  the  current  observations  shown  in  Figure  7  seem  to  imply  that  there  was  weak  flow  toward 
the  bank  beneath  the  vortex. 

The  flow  pattern  on  the  upstream  side  and  over  the  EFGB  was  also  captured  by  a  dye-release  survey  con¬ 
ducted  just  after  sections  shown  in  Figures  7a-7d  were  finished.  The  dye  patch  was  released  approximately 
at  a  depth  of  30  m  (roughly  between  23.8  and  24.1  kg  m  3  isopycnals).  It  was  sampled  by  the  ScanFish  and 
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Figure  4.  Current  velocity  components  (cm  s  '):  east/west  at  (a)  M6  and  (c)  M8  and  north/south  at  (b)  M6  and  (d)  M8;  positive  values  are  toward  east  and  north. 


the  REMUS  AUV  between  6  and  7  June  201 1  (Figure  10a).  The  dye  patch  spread  vertically  between  25  and 
45  m  depths  but  it  was  generally  trapped  within  original  density  surfaces  indicating  weak  diapycnal  mixing. 
Figure  10b  shows  a  horizontal  spreading  of  the  dye  with  the  currents.  The  currents  shown  in  Figure  10b 
(black  arrows)  are  averages  of  those  measured  between  25  and  45  m  by  shipboard  ADCPs  and  by  moored 
ADCPs  for  the  time  period  between  20:00  UTC  on  6  June  and  14:00  UTC  on  7  June  201 1.  The  dye  concentra¬ 
tion  was  normalized  by  dividing  by  0.05  ppb,  which  was  the  resolution  of  the  fluorometer.  The  current  and 
the  dye  concentration  structure  implied  flow  divergence  on  the  upstream  side  with  a  dye  patch  being 
advected  by  the  currents  around  the  southern  part  of  the  bank. 

Figures  1 1  and  12  illustrate  current  and  hydrographic  fields  when  the  mean  flow  was  westward.  Figure  1 1 
displays  currents  from  the  shipboard  ADCP  instruments  and  density  from  T8,  T9,  and  T10  transects  collected 
between  13:30  UTC  on  11  June  and  03:30  UTC  on  12  June  201 1,  while  Figure  12  shows  velocity  vectors 
along  the  transects  and  from  the  moored  ADCPs  at  four  different  depths.  There  were  similarities  in  flow 


Figure  S.  Time  series  of  (a)  temperature  (°Q,  (b)  salinity  (psu),  and  (c)  density  (kg  m  3)from  theS8  mooring;  black  dashed  lines  in  Figure 
5a  indicate  sensors  depths. 
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Figure  6.  (a)  Mean  density  (kg  m  3)  and  (b)  mean  buoyancy  frequency  (s  ')  profiles  estimated  from 
observations  recorded  at  string  moorings  (SI  58)  for  the  mean  eastward  (EF,  circles)  and  westward 
(WF,  triangles)  flow  phases. 


once  more  very  little  vertical  mixing,  at  least  along  the  ScanFish  and  REMUS  tracks 
dye  release. 


patterns  between  both 
mean  flow  phases.  On 
the  upstream  side  of  the 
bank,  the  flow  deceler¬ 
ated  and  diverged  near 
the  bank.  Similar  to  the 
eastward  flow,  the  west¬ 
ward  flow  tended  to 
generate  eddy-like  fea¬ 
tures  near  mooring  M2 
on  the  lee  side  of  the 
bank  (Figures  11c,  1 1d, 
12c,  and  1 2d).  The  sec¬ 
ond  dye  release  at  a 
depth  of  30  m  occurred 
on  8  June  201 1  on  the 
eastern  (upstream)  side 
of  the  EFGB.  The  vertical 
spreading  of  the  dye 
was  almost  negligible, 
and  the  patch  was  yet 
again  simply  adverted 
by  the  currents  over  the 
bank  just  north  of  its 
highest  part  (Figure  13). 
The  dye  experiment 
showed  that  there  was 
and  at  the  depths  of  the 


Figures  7  and  11  also  display  examples  of  potential  density  (ct8)  along  the  transects,  while  the  buoyancy  fre¬ 
quency  along  T4and  T10  is  shown  in  Figure  14.  It  is  obvious  from  the  displayed  density  and  N  distributions 
that  waters  in  the  EFGB  region  were  very  highly  stratified.  The  upper  40  m  of  the  water  column  contained  a 
number  of  well-defined  density  gradients  where  the  buoyancy  frequency  was  also  high  (N  was  often  larger 
than  0.03  s  ’).  At  depths  of  3-5  m  below  the  sea  surface,  there  was  a  first  pycnodine  (see,  for  instance, 
near-surface  density  in  Figure  7e)  partially  related  to  a  daily  temperature  cycle,  which  is  quite  pronounced 
in  the  Gulf  of  Mexico  [Stommel  et  a!.,  1 969],  and  partly  associated  with  lower-salinity  near  the  sea  surface. 
The  next  distinct  density  gradient  displayed  in  Figure  1 1  appeared  to  be  associated  with  salinity  variability 
and  was  generally  at  depths  between  10  and  15  m  depending  on  a  location.  The  main  pycnodine,  strongly 
related  to  the  temperature  gradient,  was  found  at  depths  between  35  and  40  m.  Flence,  these  data  suggest 
that  during  the  experiment  there  was  no  mixed  layer  of  significant  thickness  over  the  EFGB.  An  identical 
assessment  was  almost  true  for  the  bottom  boundary  layer  because  the  density  and  N  distributions  showed 
that  waters  were  also  stratified  near  the  bottom.  The  thickness  of  this  layer  was  estimated  using  potential 
density.  Only  waters  with  ae>  (ermax  Act)  were  considered  within  the  bottom  boundary  layer,  where  amax 
is  the  maximum  potential  density  of  each  VMP  profile  and  Act  0.01  kg  m  3.  Estimates  of  this  thickness 
from  density  data,  which  extended  nearly  to  the  bottom,  showed  that  the  bottom  boundary  layer  was  less 
than  4  m  thick  for  75%  of  the  VMP  profiles.  On  average,  this  layer  thickness  was  about  2.3  m  (standard 
deviation  1.5  m  and  median  1.4  m)  and  3.1  m  (standard  deviation  1.8  m  and  median  25m)forthe 
eastward  and  westward  mean  flow  phases,  respectively. 


5.  Turbulence  and  Mixing 

Rough  topography  on  continental  shelves  is  usually  characterized  by  enhanced  turbulence  and  mixing  [e.g., 
Loder  et  al.,  1 992;  Mourn  and  Nash,  2000;  Nash  and  Mourn,  2001 ;  Dewey  et  al„  2005].  During  the  June  experi¬ 
ment,  flow  conditions  over  the  EFGB  were  weak  to  moderate  and  stratification  was  strong.  At  the  same 
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Figure  7.  (a,  c,  and  e)  East/west  and  (h  d,  and  f)  north/south  current  velocities  (cm  s  color  bar)  from  PT5  (Figures  7a  and  7b),  PT4  (Figures  7c  and  7d),  and  T3  (Figures  7e  and  7f)  sec 
tions  taken  between  12:15  UTC  on  6  June  and  230  UTC  on  7  June  201 1  during  the  mean  eastward  flow  phase;  density  is  also  shown  (kg  m  3;  light  gray  lines  plotted  every  0.2  kg  m  3); 
top  axes  in  Figures  7a  and  7b  list  longitudes  (°;W);the  insert  shows  locations  of  the  transects  (measurements  along  PT4  and  PTSwere  taken  from  the  R/V  Pelican;  T3  data  were  collected 
from  the  IW  Manta). 

time,  however,  observed  turbulence  and  mixing  over  the  bank  and  on  its  lee  side  were  elevated.  Figures  15 
and  1 6  display  turbulent  kinetic  energy  dissipation  rates  («)  for  the  eastward  and  westward  mean  flows  for 

selected  transects,  respectively.  The  dissipation  rates  were  estimated  from  e  v(%-)2  assuming  that  tur¬ 
bulence  was  isotropic.  This  assumption  holds  when  the  ratio  of  e/vN2  is  larger  than  200  [Gargett  et  al.,  1984]. 
In  case  of  dissipation  rates  presented  in  Figures  15and  16,  the  ratio  was  larger  than  200  for  14.3%  of  the  e. 
estimates.  Itsweire  et  al.  [1993]  concluded  that  an  assumption  of  isotropy  may  lead  to  underestimation  of  e 
by  factors  2-4  when  the  ratio  is  less  than  200  and  the  Richardson  number  (Ri  N 2  (^)  2,  where  ^  is  the 
horizontal  current  shear)  is  larger  than  0.37  which  was  the  case  during  the  experiment  (more  discussion 
concerning  the  Richardson  numbers  estimated  from  the  MORT  observations  can  be  found  at  the  end  of  this 
section).  That  means  that  for  57.6%  of  the  dissipation  rate  estimates  are  underestimated.  Itsweire  etal. 

[1 993]  also  suggested  that  turbulence  produces  no  significant  net  buoyancy  flux  and  insignificant  mixing  if 
the  same  ratio  is  less  than  19,  which  is  the  case  for  28.1%  of  the  estimates  here. 

The  upstream  side  of  the  bank  was  better  sampled  during  the  eastward  mean  flow  phase  by  the  microstruc¬ 
ture  measurements  (Figure  15).  Except  very  near  the  sea  surface  or  the  bottom,  e  was  generally  less  than 
5X10  8  W  kg  1  on  the  upstream  side  with  depth-averaged  means  and  medians  off:  below  10  8  W  kg  1 
for  each  profile.  Away  from  the  sea  surface  and  bottom,  spatial  structures  of  the  dissipations  rates  also  indi¬ 
cate  that  turbulence  was  rather  intermittent  in  the  water  column  over  the  bank  and  on  its  downstream 
side.  Some  of  those  turbulent  patches  on  the  lee  side  may  have  been  assoaated  with  possible  vortices  as,  for 
instance,  the  patches  found  between  depths  of  40  and  80  m  starting  at  about  3  km  on  the  lee  side  of  T3  tran¬ 
sect  where  e  was  as  high  as  10  6Wkg  '  (Figure  15b).  High  values  of  e,  which  easily  reached  10  sWkg 
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Figure  8.  Current  velocities  (cm  s  1 )  dlong  PT5,  PT4,  and  T3  at  (a)  20  m,  (b)  52  m,  (c)  60  m,  and  (d)  76  m  as  measured  by  shipboard  (blue 
arrows)  and  moored  (cyan  arrows)  ADCP  instalments;  also  shown  are  time  velocity  means  (12:15  UTC6  June  to  2:30  UTC  7  June  2011; 
green  arrows)  estimated  from  mooring  observations,  and  depth  contours  for  the  chosen  depth  levels  are  highlighted  by  black  contour 
lines. 


were  also  found  near  the  bottom  of  the  EFGB  as  well  as  near  the  sea  surface  (above  10  m;  Figures  15  and 
1 6).  Near-surface  turbulence  seemed  to  be  partially  related  to  the  daily  temperature  cycle  since  winds 
were  rather  weak  and  surface  waves  were  small  and  not  breaking  (surface  wave  significant  heights  were 
less  than  0.5  m  [Wijesekera  et  at.,  2013]).  Additionally,  estimates  of  the  Kolmogorov  scale  (L0  e,/2N  3/2) 
indicated  that  length  scales  of  turbulent  overturns  in  the  EFGB  region  were  less  than  1  m  (99%  of  L0 
estimates  <  1 ). 

Measured  rates  of  energy  dissipation  and  mixing  by  turbulence  per  unit  width  over  the  bank  and  on  its 
downstream  side  were  calculated  from  J  J  p(e,+Jb)dzdx,  where  Jb  is  the  irreversible  buoyancy  flux  that 
is  often  approximated  by  Te  (f  is  the  mixing  efficiency  and  it  is  often  assumed  to  be  0.2  [Osbom,  1980]). 
They  varied  highly  among  transects.  Integrations  were  done  vertically  between  the  first  and  last  depth 
where  energy  dissipation  rates  were  available  and  then  horizontally  between  the  foot  of  the  bank  on  the 
upstream  side  to  the  last  available  VMP  profile  on  the  lee  side  of  the  EFGB.  For  those  taken  during  the 
mean  eastward  flow  and  shown  in  Figure  15,  they  were  28.1, 28.5,  and  65Wm  '  for  T2,  T3,  and  T4,  respec¬ 
tively.  Estimates  for  the  T8,  T9,  and  T10  transects,  displayed  in  Figure  16,  were  correspondingly  13.2, 40.0, 
and  23.8  W  m  ’.  They  also  varied  temporarily  at  the  same  location  and  for  the  same  mean  flow  conditions. 
For  instance,  for  the  most  northward  transect  (indicated  as  T10  in  Figure  16),  they  were  23.8  and  9.7  W  m  ' 
as  estimated  from  observations  taken  about  10  h  apart  on  12  June  2011. 

Following  Osbom  [1980],  the  eddy  diffusivity  (Kp)  was  estimated  from  Kp  r i:N  2.  It  has  been  shown  that 
the  mixing  efficiency,  T,  may  not  be  completely  constant  particularly  when  stratification  is  weak  [e.g,  Shih 
et  al„  2005;  Ivey  et  at.,  2008;  Mashayek  et  al.,  201 3;  Bluteau  et  at.,  20 1 3],  For  eddy  diffusivity  estimates  dis¬ 
cussed  here,  T  0.2  was  used  as  presented  in  the  vast  majority  of  oceanographic  mixing  literature  [see,  for 
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Figure  9.  (a)  Temperature  ( V;°Q  and  (b)  density  (<?„';  kg  m  3)  anomalies  along  PT4,  (c)  relative  vorticity  (s  ')  estimated  from  velocity 
observations  from  PT4  and  PT5;  H  is  the  water  depth. 


example,  Osborn,  1980;  Oakey,  1982 ;  Mourn,  1996;  Nash  and  Mourn,  2001].  Gregg  et  al.  [2012]  also  provide  an 
extensive  discussion  and  justification  for  this  choice  of  T.  Additionally,  the  eddy  diffusivity  was  estimated 
only  for  N2  >  1 0  6  s  2  to  avoid  a  bias  toward  large  values  of  K(,  at  depths  where  stratification  was  very 
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Figure  10.  (a)  The  location  (red  square)  of  the  dye  release  on  6  June  2011  and  the  REMUS  (blue)  and  ScanFish  (magenta)  survey  tracks 
from  6  to  7  June  201 1;  (b)  normalized  dye  concentration  (color  bar)  and  currents  (cm  s  black  arrows). 
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Figure  11.  (a,c,  and  e)  East/west  and  (b,  d,  and  f)  north/south  current  velocities  (cm  s  '.color  bar)  fromTIO  (Figures  1  la  and  11b),  T9  (Figures  11c  and  1 1d),  and  T8  (Figures  lie  and 
1 1f)  sections  taken  between  13:45  UTC  on  1 1  June  and  3:30  UTC  on  12  June  201 1  during  the  mean  westward  flow  phase;  density  is  also  shown  (kg  m  3;  light  gray  lines  plotted  eveiy 
0.2  kg  m  3);top  axes  in  Figures  11a  and  1 1b  list  longitudes  C  ;  W);  the  insert  shows  locations  of  the  transects  (measurements  along  all  transects  were  collected  from  the  R/V  Manta). 

weak.  Estimates  of  Kp  spanned  from  2.2  X  10  7  to  1.1  X  10  2  m2  s  '  and  varied  horizontally  and  verti¬ 
cally.  The  horizontal  variations  of  the  eddy  diffusivity  along  with  means  and  medians  for  individual  pro¬ 
files  are  shown  in  Figure  1 7  for  the  T3  and  T9  transects.  Note  that  when  the  approach  proposed  by  Shih 
et  al.  [2005]  is  employed  and  e/vN2>  1 00  (23.1%  cases),  estimates  of  the  eddy  diffusivity  are  up  to  an 
order  of  magnitude  lower  than  those  shown  in  Figure  17.  Upstream  of  the  EFGB,  the  mixing  was  rather 
subdued  with  Kp  generally  less  than  5X10  5  m2  s  \  Such  a  low  values  of  the  eddy  diffusivity  are  repre¬ 
sentative  for  mixing  over  slowly  varying  topography  in  the  ocean  and  in  the  open  ocean  thermocline. 
Higher  Kp  (>10  4  m2  s  ’)  and  more  vigorous  mixing  was  found  over  and  downstream  of  the  bank.  Over 
the  EFGB,  estimates  of  the  eddy  diffusivity  were  10  3  m2  s  1  or  higher,  and  elevated  mixing  was  often 
found  near  the  bottom.  On  the  lee  side,  similar  high  Kp  values  were  found  at  depths  where  the  velocity 
observations  implied  a  presence  of  vortices.  Spatial  variability  of  mixing  in  the  EFGB  region  was  also 
reflected  by  results  of  the  dye  experiment.  For  both  releases,  the  dye  patches  barely  spread  vertically 
implying  weak  vertical  mixing  as  indicated  by  rather  low  estimates  of  the  eddy  diffusivity,  especially  on 
the  upstream  side  of  the  EFGB.  The  dye  patches  were  advected  horizontally  by  the  currents  generally  at 
depths  where  they  were  injected. 

Despite  spatial  variations  of  the  eddy  diffusivity,  turbulent  mixing  was  still  more  effective  and  stronger  over 
the  EFGB  when  compared  to  that  over  the  flat  or  slowly  varying  bottom  of  the  continental  shelf.  One  way 
to  compare  mixing  proficiency  is  to  evaluate  total  diapycnal  buoyancy  fluxes  over  the  bank  and  shelf  areas, 
i.e.,  JtAb  and  JSAS  where  Jb,  Js  are  the  fluxes  and  Ab,  As  are  the  areas  of  the  bank  and  shelf,  respectively.  The 
fluxes  are  defined  as  Jb  KflbN*  and  Js  KpsN2,  where  K^,  Kps  are  eddy  diffusivities  and  Wjj,  N2  are  squared 
buoyancy  frequencies  for  the  bank  and  the  shelf,  respectively.  Following  Lueck  and  Madge  [1 997],  the  area 


JAROSZETAL. 


©2014.  American  Geophysical  Union.  All  Rights  Reserved. 


6414 


AGU  Journal  of  Geophysical  Research:  Oceans 


10.1 002/201 4JC009934 


Figure  12.  Current  velocities  (cm  s  ')  along  T8,  T9,  and  TIOat  (a)  16  m,  (b)  32  m,  (c)  52  m,  and  (d)  64  mas  measured  by  shipboard  (blue 
arrows)  and  moored  (cyan  arrows)  ADCP  instruments;  also  shown  are  time  velocity  means  (13:45  UTC 1 1  June  to  3:30  UTC  12  June  2011; 
green  arrows)  estimated  from  mooring  observations,  and  depth  contours  for  the  chosen  depth  levels  are  highlighted  by  black  contour 
lines. 

of  a  flat-bottom  shelf,  which  is  required  to  perform  as  efficient  mixing  as  that  over  the  EFGB,  can  be  esti¬ 
mated  by  setting  JbAb  J^s.  The  area  of  the  EFGB  within  a  100  m  isobath  is  about  40.76  km2.  Taking  a  typi¬ 
cal  value  of  the  eddy  diffusivity  for  the  shelf,  Kps  10  Sm2s  and  those  estimated  for  the  EFGB, 
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Figure  13.  (a)  The  location  (red  square)  of  the  dye  release  on  8  June  2011  and  the  REMUS  (blue)  and  ScanFish  (magenta)  survey  tracks 
from  8  to  9  June  201 1;  (b)  normalized  dye  concentration  (color  bar)  and  currents  (cm  s  black  arrows). 


JAROSZETAL. 


©2014.  American  Geophysical  Union.  All  Rights  Reserved. 


6415 


AGU  Journal  of  Geophysical  Research:  Oceans  10.1002/2014JC009934 


-93.62 


-93.6 


-93.58 


f 

I 


0.04 

0.03 

0.02 

0.01 

0 


/b)  - 1 - 


20 
1"  40 


■5.  60 

a> 

°  80 


100 


-1  0 
Distance  (km) 


Figure  14.  Buoyancy  frequency  (s  ')  along  (a)  T4and  (b)  T10  transects;  the  insert  shows  locations  of  the  transects. 
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Figure  IS.  Turbulent  kinetic  energy  (TKE)  dissipation  rates  (j;  W  kg  ')  along  (a)T4,  (b)  T3,  and  (c)  T2  transects;  the  top  axis  in  Figure  15a 
lists  longitudes  f;  W);  the  insert  shows  locations  of  the  transects;  microstructure  observations  were  taken  between  21  f)0  UTC  on  5  June 
and  00:30  UTC  on  8  June  201 1  from  the  R/V  Manta  during  the  mean  eastward  flow  phase;  locations  of  the  VMP  drops  are  marked  by  black 
dots. 
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Figure  16.  Turbulent  kinetic  energy  fTKE)  dissipation  rates  (e;  W  kg  ')  along  (a)T10,  (b)  T9,  and  (c)T8  transects;  the  top  axis  in  Figures  16a 
lists  longitudes  p;  W);  the  insert  shows  locations  of  the  transects;  microstructure  observations  were  taken  between  13:45  UTCon  11  June 
and  3:30  UTC  on  1 2  June  201 1  from  the  R/V  Manta  during  the  mean  westward  flow  phase;  locations  of  the  VMP  drops  are  marked  by  black 
drops. 


10  4  1 0  3  m2  s  \  and  assuming  that  Nl  ~  N 2  the  shelf  area  has  to  be  from  10  to  100  times  larger 
than  that  of  the  bank,  i.e.,  it  should  be  between  407.6  and  4076  km2. 

The  velocity  and  hydrographic  observations  taken  by  the  shipboard  ADCP  and  the  VMP  profiler  were  also 
used  to  calculate  the  Richardson  number  (estimates  not  shown)  along  transects  shown  in  Figures  15  and 
1 6.  Values  of  Ri  were  rather  underestimated  due  to  coarse  vertical  resolution  of  the  current  velocity  observa¬ 
tions  (4  m);  thus,  higher-resolution  current  observations  are  required  to  obtain  a  more  accurate  spatial  Ri 
distribution  and  to  fully  understand  mixing  over  the  EFGB.  Additionally,  the  Richardson  number  was  esti¬ 
mated  only  for  the  squared  current  shear  and  N2  larger  than  10  6s  2  to  avoid  a  bias  toward  very  small  val¬ 
ues  of  Ri  at  depths  where  the  current  shear  was  low  or/and  stratification  was  weak.  The  distribution  of  Ri 
does  not  indicate  that  mixing  in  the  EFGB  region  was  mainly  related  to  shear  instabilities,  i.e.,  the  Richard¬ 
son  number  was  rarely  found  to  be  0.25  or  less.  Ri  was  generally  above  1  on  the  upstream  side.  The  smaller 
values  of  the  Richardson  number  at  some  depths,  even  as  low  as  0.25  or  less,  were  mostly  found  on  the 
downstream  side  and  over  the  crest  of  the  bank. 

6.  Discussion 

Flows  over  obstacles  are  often  described  in  terms  of  dimensionless  parameters.  Among  those  parameters, 
more  common  ones  are:  the  Froude  number  (Fa  nU/NH),  a  quantity  often  referred  to  as  the  topographic 
Froude  number  {Fh  U/Nh),  the  Rossby  number  [Ro  U/fL),  the  nondimensional  obstacle  width  (/„  Nl/U) 
and  height  (hn  h/H)  where  h  is  the  obstacle  height,  /  is  the  half  of  the  obstacle  width  (L),  H  is  the  finite 
upstream  depth,  U  is  the  upstream  flow  speed,  and  f  is  the  Coriolis  parameter  [see,  for  example,  Hunt  and 
Snyder,  1 980;  Baines,  1 995,  Vosper  et  at.,  1 999],  Moreover,  the  ratio  (an)  of  the  flow  depth  (H  1 00  m)  to  the 
obstacle  width  (L  6000  m)  is  a  good  measure  whether  the  flow  is  hydrostatic.  If  an  is  much  less  than  one, 
i.e.,  for  an  obstacle  with  a  long  horizontal  scale,  the  flow  is  assumed  to  be  hydrostatic,  which  is  a  good 
approximation  for  the  EFGB  for  which  this  ratio  is  ~  10  2.The  rotational  effects  were  also  important  in  this 
region  since  Ro  was  033  (the  Ro  standard  error  0.06;  f  6.84X  10  5s  ’). 
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Figure  17.  (a  and  b)  Eddy  diffusivities  (m2  s  ')  and  (c  and  d)  bathymetry  for  T3  (Figures  17a  and  17c)  and  19  (Figures  17b  and  17d)  trans 
ects;  means  and  medians  of  eddy  diffusivities  for  each  profile  are  marked  by  red  and  blue  lines,  respectively. 


The  nondimensional  obstacle  width,  which  is  larger  than  100  for  the  EFGB,  determines  ifthe  flow  produces 
dispersive,  nonhydrostatic  waves  (/„  <  10)  or  hydrostatic  waves  confined  mostly  to  the  region  near  the 
obstacle  (f„  >  1 0).  This  number  indicates  that  the  latter  scenario  was  a  possibility  in  the  EFGB  region,  and 
the  density  structure  (e.g.  Figure  7e)  seems  to  indicate  that  wave-like  features  were  generated  over  the 
bank  top.  Additionally,  the  nondimensional  height  (fin)  of  the  EFGB  varies  but  is  not  larger  than  0.82  (a  h 
maximum  ~  82  m)  and  for  a  major  part  of  the  bank,  i.e.,  for  82%  of  the  EFGB  area,  h„  is  below  0.5.  hn  is  not 
likely  to  be  important  unless  it  is  very  close  to  1.  In  such  a  case,  the  free  surface  may  interact  with  an  obsta¬ 
cle  and  generated  disturbances  are  confined  vertically.  For  the  EFGB,  that  may  be  the  case  only  at  its  very 
small  southern  part  where  the  bank  is  over  80  m  high.  Previous  research  of  the  flow  over  a  three- 
dimensional  obstacle  also  has  shown  that  the  response,  especially  on  a  lee  side  of  the  obstacle,  strongly 
depends  on  its  shape  and  asymmetry  [e.g,  Baines,  1995],  The  observations  presented  here  imply  that  inter¬ 
actions  between  the  free-stream  flow  and  the  asymmetric  EFGB  were  very  complex  during  the  field 
experiment 

The  Froude  number  can  be  used  to  classify  flows  over  topography.  In  ocean  science,  the  classically  defined 
Froude  number  is  based  on  the  first-mode  internal  wave  speed  (cj  (e.g,  Fa  U/ca  nU/NH,  where  ea 
NH/n).  If  F„  >  1,  linear  internal  waves  do  not  propagate  upstream,  while  if  Fa  <  1  it  is  possible  that  the  first 
wave  mode  (at  least)  may  propagate  upstream  as  a  shear  front  or  columnar  disturbance  mode  [e.g,  Wei 
et  a/,  1 975;  Baines,  1977].  It  is  also  utilized  for  flow  descriptions  in  the  hydraulic  sense.  Lawrence  [1 993]  has 
distinguished  four  basic  flow  regimes:  subcritical,  crest-controlled,  approach-controlled,  and  supercritical. 
The  flow  is  referred  to  as  subcritical  if  F0  is  very  low  (<<1).  Transitional  crest-controlled  and  approach- 
controlled  flows  occur  in  intermediate  Fa  (~1),  while  supercritical  flows  are  characterized  by  F„>  1.  During 
the  time  period  (4-12  June  201 1)  of  the  extensive  sampling  over  the  EFGB,  hourly  estimates  of  Fa  were 
always  below  051 .  For  the  eastward  mean  flow,  its  average  values  (averaged  over  a  time  used  to  collect 
observations  along  a  transect)  and  their  standard  errors  (in  parentheses)  were  0.38  (±0.02),  0.32  (±0.01 ), 
and  0.25  (±0.01)  forT2,  T3,  andT4,  respectively.  All  Fa  values  were  below  0.20 for  the  westward  mean  flow, 
e.g,  the  transect  means  were  0.11  (±0.01),  0.09  (±0.01),  and  0.17  (±0.02)  forT8,T9,  and  T10,  respectively.  H 
was  taken  as  100  m  and  estimated  40  m  depth  averages  of  N  and  U  were  used  in  the  Fa  calculations.  These 
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transect-averaged  Froude  numbers  were  much  less  than  1  indicating  that  the  flow  was  subcritical  during 
the  experiment  in  the  EFGB  region. 

Estimated  mean  topographic  Froude  numbers,  which  are  commonly  used  in  atmospheric  research  to  study 
airflows  over  mountainous  terrains,  were  less  than  0.25  for  all  transects  regardless  of  the  mean  flow  phase, 
and  these  small  numbers  clearly  imply  that  the  flow  was  also  subcritical  and  highly  nonlinear.  It  has  been 
shown  that  highly  stratified  flows  with  small  Fa  and  Fh  numbers  are  spatially  and  temporally  variable  with 
nonlinear  effects  being  important  [e.g,  Baines,  1979;  Hunt  and  Snyder,  1 980;  Smolarkiewicz  and  Rotunno, 

1 989;  Ding  et  at.,  2003]  when  compared,  for  instance,  to  weakly  stratified  flows  with  Froude  numbers  larger 
than  1  [e.g.,  Smith,  1988].  Moreover,  these  subcritical-highly  stratified  flows  are  blocked  on  the  upstream 
side  of  the  obstacle  and  tend  to  move  around  it  [e.g,  Brighton,  1 978;  Baines,  1979;  Hunt  and  Snyder,  1980]. 
This  is  the  case  here,  i.e,  even  though  the  observations  were  spatially  very  coarse  they  still  were  able  to  cap¬ 
ture  the  decelerated  and  diverged  flow  on  the  upstream  side  of  the  EFGB  (Figures  7, 8, 11,  and  12).  Shep¬ 
pard  [1956]  postulated  that  fluid  parcels  originated  at  depth  levels  above  that  of  h(l  -  Fh)  would  pass  over 
the  top.  Unfortunately,  the  observations  did  not  resolve  exact  depth  levels  where  a  separation  between 
fluid  parcels  going  over  and  those  flowing  around  occurred  for  the  entire  EFGB  during  the  experiment; 
however,  applying  this  approach  to  T3  and  PT4  transects,  only  fluid  at  depths  above  0.8fi  (~44  m)  would 
pass  over  the  bank,  while  fluid  at  and  below  that  depth  would  move  around  the  bank  (the  bank  height  (h) 

~  70  m).  The  current  velocity  observations  from  these  two  transects  show  that  the  flow  was  over  the  bank 
approximately  to  40  m  below  the  sea  surface,  i.e„  just  above  the  depth  estimated  from  Sheppard's  criterion, 
whereas  the  flow  decelerated  and  began  moving  along  the  bank  below  40  m. 

Previous  laboratory  and  numerical  studies  [e.g,  Baines,  1 979;  Hunt  and  Snyder,  1 980;  Smolarkiewicz  and  Rotunno, 
1 989;  Schar  and  Durran,  1 997;  Vosper  et  al„  1 999]  have  shown  that  a  wake  filled  with  vortices  develops  on  the 
lee  side  of  an  obstacle  for  subcritical  flows.  The  lee  wake  is  usually  characterized  by  elevated  turbulence  levels 
and  mixing  when  compared  to  those  observed  on  the  upstream  side  of  the  obstacle.  The  studies  have  also 
shown  that  vortex  shedding  occurs  generally  when  Fh  falls  below  0.4  [Vosper  et  al„  1999],  and  that  shed  vortices 
are  considered  to  be  a  means  of  dissipating  flow  energy  [Schar  and  Durran,  1 997].  Our  current  veloaty  measure¬ 
ments  and  estimated  relative  vorticity  indicate  that  there  was  a  cyclonic  vortex  on  the  lee  side  (the  near  M4  and 
M4a  moorings;  see  Figure  1  for  their  locations)  of  the  EFGB  during  the  mean  eastward  flow.  Results  from  numer¬ 
ical  simulations  by  Schar  and  Durran  [1 997]  show  that  all  shed  cydonic  and  anticydonic  vortices  are  warm-core 
vortices,  which  may  explain  positive  temperature  anomalies  observed  in  the  cydonic  vortex  near  M4  and  M4a. 
An  eddy-like  structure  emerging  from  the  current  observations  was  less  clear  for  the  mean  westward  flow;  how¬ 
ever,  the  data  seem  to  imply  that  a  vortex  was  generated  on  the  lee  side  of  the  bank  near  M2  (see  Figure  1  for 
the  mooring  location).  Our  observations  are  not  sufficient  to  speculate  whether  vortices  observed  during  the 
subcritical  flow  conditions  in  the  EFGB  area  were  consistent  features  that  were  generated  by  the  mean  currents, 
or  if  they  were  periodic  features  forced  by  tidal/inertial  flows,  or  if  they  resulted  from  a  combination  of  both 
motions.  Moreover,  the  collected  microstructure  data  dearly  show  that  turbulence  and  mixing,  as  expected, 
were  enhanced  on  the  downstream  side  when  compared  to  those  on  the  upstream  side  of  the  EFGB  for  both 
mean  flow  phases  (Figures  15-17).  Enhanced  turbulence  was  patchy  and  was  partly  related  to  vortices. 

Furthermore,  recent  analyses  of  the  pressure  distribution  by  Wijesekera  etal.  [2014]  have  shown  that  the  highly 
stratified,  nonlinear,  hydrostatic  flow  generated  high  form  drag  over  the  EFGB.The  form  drag  occurs  when  cur¬ 
rents  flow  over  rough  topography  creating  a  pressure  difference  between  the  upstream  and  downstream  sides  of 
an  obstacle,  i.e,  a  force  that  opposes  the  flow  [e.g.  Mourn  and  Nash,  2000;  Edwards  et  al.,  2004;  Warner  et  al.,  201 3; 
Warner  and  MacC ready,  201 4],  This  drag  ranged  from  about  2000-3000  N  m  1  and  shows  multiple  time  scale  vari¬ 
ability  over  the  EFGB.  Similarly,  estimated  average  frictional  drag  over  the  bank  was  also  quite  high  (~0.006). 
These  findings  are  in  agreement  with  those  reported  for  other  underwater  banks  [e.g.  Mourn  and  Nash,  2000]. 

7.  Conclusions  and  Summary 

Beginning  on  4  June  201 1,  a  9  day  extensive  oceanographic  survey  was  conducted  over  the  East  Flower 
Garden  Bank,  a  coral  reef,  located  in  the  northwestern  Gulf  of  Mexico  about  190  km  southeast  of  Galveston, 
TX.  The  EFGB  is  considered  a  rough  topographic  feature.  It  is  about  6  km  wide,  10  km  long,  and  is  located  at 
the  shelf  edge  in  100  m  of  water  depth  with  a  peak  rising  to  about  18  m  below  the  sea  surface.  The  field 
experiment  involved  two  research  vessels,  the  R/V  Manta  and  the  R/V  Pelican.  Observations  were  collected 
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by  shipboard  (ADCP),  ship-towed  (ScanFish),  ship-deployed  (VMP),  autonomous  (AUV),  and  moored  (TCP 
string  and  ADCP  moorings)  instrumentation.  These  data  were  then  used  to  characterize  the  flow  variability, 
stratification,  turbulence,  and  mixing  over  the  surveyed  bank. 

During  the  experiment  the  mean  currents  over  the  EFGB  switched  from  eastward  to  westward  on  8  June 
201 1.  The  measured  flow  was  a  combination  of  the  subtidal,  inertial,  tidal,  and  higher-frequency  motions 
with  a  maximum  speed  of  about  50  cm  s  ’.  The  waters  were  highly  stratified  throughout  the  entire  water 
column  with  the  depth  mean  buoyancy  frequency  of  0.016  and  0.018  s  1  for  the  eastward  and  westward 
mean  flow  phases,  respectively.  As  a  result  the  mixed  surface  and  bottom  boundary  layers  were  generally 
less  than  4  m  thick.  Consequently,  the  flow  conditions  over  the  EFGB  were  classified  as  subcritical  when 
transects  were  sampled  since  all  estimated  Froude  numbers  (F0)  were  less  than  0.4.  Based  on  computed 
nondimensional  parameters  ( an ,  /„,  Fh,  hn,  and  (to),  flow  dynamics  over  the  bank  were  also  hydrostatic,  non¬ 
linear,  and  rotational  effects  were  important. 

Observations  clearly  showed  that  flow  structure  and  mixing  were  highly  dependent  on  the  direction  and 
strength  of  the  currents  in  the  EFGB  region;  thus,  they  varied  spatially  and  temporarily.  As  past  observatio¬ 
nal,  laboratory,  and  numerical  efforts  [e.g.,  Baines,  1 979;  Hunt  and  Snyder,  1 980;  Smolarkiewicz  and  Rotunno, 

1 989;  Vosper,  2000;  Nash  and  Mourn,  2001;  Skyllingstad  and  Wijesekera,  2004;  Dewey  et  al.,  2005]  have  indi¬ 
cated,  responses  resulting  from  interactions  between  the  free-stream  flow  and  the  bank  were  significantly 
different  on  the  upstream  and  lee  sides  of  the  EFGB.  For  highly  stratified,  subcritical  flows  over  an  obstacle, 
blocking  and  diverging  of  the  flow  is  expected  on  its  upstream  side,  and  these  flow  features  were  observed 
on  the  upstream  side  of  the  EFGB  for  both  mean  eastward  and  westward  current  phases.  On  the  down¬ 
stream  side  of  the  bank,  a  wake  developed  as  identified  by  past  studies  of  flows  with  low  Froude  numbers. 
The  current  observations  and  estimated  relative  vorticity  also  implied  that  eddy-like  features  were  present 
on  the  lee  side  of  the  bank.  Furthermore,  turbulence  was  amplified  over  the  EFGB  top  and  on  its  lee  side. 
Estimated  TKE  dissipation  rates  from  small-scale  velocity  shear  observations  were  as  high  as  10  6  W  kg  1 
resulting  in  measured  rates  of  energy  dissipation  and  mixing  by  turbulence  per  unit  width  over  the  bank  as 
high  as  40  W  m  \  Even  for  such  weak  flow  conditions,  mixing  on  the  downstream  side  was  also  enhanced 
by  a  few  orders  of  magnitude  above  the  typical  shelf  value  of  10  5  m2  s  \  and  the  eddy  diffusivity  there 
was  as  high  as  10  3  m2  s  ’.On  the  upstream  side,  estimated  values  of  KfJ  were  close  to  that  for  the  conti¬ 
nental  shelf  with  gradually  varying  bathymetry  and  were  generally  less  than  0.5  X  10  4  m2  s  \ 

Significantly  stronger  flows  over  the  EFGB  were  recorded  by  the  year-long  moorings  (Ml  -M5)  [ Teague  et  al., 
2013].  Turbulence  and  mixing  are  expected  to  be  more  elevated  during  the  more  extreme  currents  com¬ 
monly  observed  in  the  EFGB  region  than  turbulence  and  mixing  observed  during  our  9  day  intense  oceano¬ 
graphic  survey.  Flence,  the  coral  reefs  can  be  considered  as  "hot  spots"  for  mixing  on  the  shelf  in  the 
northwestern  Gulf  of  Mexico. 
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